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Abstract 

We consider the problem of coloring a grid using k colors with the restriction that in each 
row and each column has an specific number of cells of each color. In an already classical 
result, Ryser obtained a necessary and sufficient condition for the existence of such a coloring 
when two colors are considered. This characterization yields a linear time algorithm for 
constructing such a coloring when it exists. Gardner et al. showed that for k > 7 the problem 
is NP-hard. Afterward Chrobak and Diirr improved this result, by proving that it remains 
NP-hard for k > 4. We solve the gap by showing that for 3 colors the problem is already 
NP-hard. Besides we also give some results on tiling tomography problems. 

1 Introduction 

Tomography consists of reconstructing spatial objects from lower dimensional projections, and has 
medical applications as well as non-destructive quality control. In the discrete variant, the objects 
to be reconstructed are discrete, as for example atoms in a cryst aline structure, see [Tj. 

One of the first studied problem in discrete tomography involves the coloring of a grid, with 
a fixed number of colors with the requirement that each row and each column has a specific total 
number of entries of each color. 

More formally we are given a set of colors C, and an mxn matrix M, whose items are elements 
of C. The projection of M is a sequence of vectors r'^ e N™, s'^ E N", for c e C, where 

r'r = \{j:M,,=c}\, ,s^= = |{»:Af,,=c}|. 

In the reconstruction problem, we are given only a sequence of vectors satisfying for 1 < i < m, 
1 < J < n, c e C, 

C C J j 

and the goal is to compute a matrix M that has the given projections. If there are fc = |C| colors, 
we call it the /c-color Tomography Problem. 

It was known since long time, that for 2 colors, the problem can be solved in polynomial 
time [8 . Ten years ago it was shown that the problem is NP-hard for 7 colors [5 . By NP- 
hardness, we mean that the decision variant — deciding whether a given instance is feasible, i.e. 
admits a solution — is NP-hard. Shortly after this proof was improved to show NP-hardness for 
4 colors, leaving open the case when \C\ = 3 [5]. This paper closes the gap, by showing that for 3 
colors already the problem is NP-hard. 
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Just to fix the notation, for \C\ 2 we denote the colors as hlack and white, and use symbols 
B,W. For \C\ = 3 we denote the colors as red, green and yellow and use symbols R,G,Y. Notice 
that we can think white and yellow as ground colors in the 2 and 3— color problem, respectively. 
Thus when we denote the instance of the tomography problem, we sometimes omit the white or 
yellow projections as they are redundant. In addition for a 2-color instance (r^,s^) we omit the 
superscript when the context permits it. 

First we recall some well known facts about the 2-COLOR tomography problem. 

Lemma 1 ([8j) Let (r, s) be a feasible instance of the 2-color tomography problem. Let I be some 
set of rows, and J be some set of columns. Lf 

^r,-^., = |/x J|, (2) 

is/ 

then every solution to the instance will be all black in I x J and all white in I x J . 

Proof: The sets I, J divide the grid into four parts, IxJ, IxJ, IxJ and I x J. The value 
X^ie/ '^^ equals the number of black cells in the first two parts, and the number of black 

cells in the second and last part. So the difference is the number of black cells in / x J minus the 
number of black cells in / x J. So when ^ holds, the first part must be all black and the last 
part all white. □ 

Before stating the next lemma, we need to introduce some notation about vectors. The 
conjugate of a vector s E {0,1,..., m}" is defined as the vector s* E {0,1,..., n}™ where 
■^i ~ \{j '■ ^ There is a very simple graphical interpretation of this. Let be an m x n 
matrix M , such in column j, the first Sj cells are colored black and the others are colored white. 
Then the conjugate of s is just the row projection of M, see figure [T] 
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Figure 1: Example of a vector s and its conjugate s* . 



Note that s* is always a non-increasing vector. If in addition s is non-increasing we have that 
(s*)* = s since in this case s* = max{j : Sj > i} and s* > j if and only if Sj > i. 

For every s,t G N" we say that s dominates t, denoted s ^ i, if X]j=i — ^^^^ every 

1 < i < n. For any < fc < n we define the set Xn,k '■— {x E {0, 1}" : — Clearly >z 

defines a partial order on Xn,k, and we show now that it has a small depth. 

Lemma 2 ([2j) Let n,k be two integers with < k < n. Suppose we have a strictly increasing 
sequence 

b^ <b^ < ...<b\ 
of vectors from Xn,k- Then q < k{n — k) . 

Proof: For each vector a E X^ ^ we associate the number (p{a) defined by f{a) — Si=i '^i- 

If a ~< P then X]^=i '^j — Pi f*-"^ every 1 < 1 <n and the inequality is strict for at least 

one t. We conclude that a < (3 implies V'(q^) < 

Therefore the vectors with extreme values for if are a — (0, . . . , 0, 1, . . . , 1) and (3 — (1, . . . , 1, 0, . . . , 0). 

fe k 

Since 1^9(0;) — k{k — l)/2 and f{/3) — k{k — l)/2 + k{n — k), this concludes the proof. □ 

A well-known characterization of the feasible instances of the 2-COLOR tomography prob- 
lem can be expressed using dominance. 
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Lemma 3 ([§]) Let (r, s) be an instance of the 2-color tomography problem, such that r is non- 
increasing. Then (r, s) is feasible if and only if r < s* . Moreover if r = s* , then there is a single 
solution, namely the realization having the first Sj cells of column j colored black, and the others 
white. 

There is a very simple graphical interpretation of this. Again let M be a matrix where in column 
j the first Sj cells are colored black and the remaining cells white. Then the row projection of M 
is s* , and if s* = r we are done. Now if s* ^ r, then some of the black cells in M have to be 
exchanged with some white cells in the same column but a lower row. These operations transform 
the matrix in such a way, that the new row projection is dominated by s* . So if s* does not 
dominate r, then there is no solution to the instance. 



2 The gadget 



The gadget depends on some integers n,k,u,v with 1 < k,u,v < n and u ^ v as well as on two 
vectors a, P G Xn,k- It is defined as the instance of n rows, and 2n + 2 columns with the following 
projections for 1 <i,j <n 



1 



if i G {u, v} 
otherwise 



if i G {u, v} 
1 otherwise 



'n+l 



°7l + 2 



n — j + aj 
1 

71 - fc + 1 





G 



'n+2 





n- 1 
fc- 1 
n — j + I 



Lemma 4 // the instance above is feasible then a ^ (3. Moreover, if a = (3 then the instance is 
feasible if and only if + > 1 . 

Proof: Assume the instance is feasible, we will show that this implies a ^ (3. Consider the yellow 
projection vectors r^ = 2n ~\- 2 — r^ — r'^ and s^ — n — s^ — s^. We have that rf = 2{n — i) + 1 
for 1 < i < n. Note that r^ is a non-increasing vector. Similarly, we obtain that sj=j — aj and 
^n+2+j — j ~ ^ + 0j ^ for 1 !i J !i ^1 and s^_^_i = s^_|_2 — 0. The conjugate of the column yellow 
projections is a vector (s^)* with 

(s^)* ^ 2{n - i) + 1 - a, + 13,. 

Then clearly ^ (s^)* if and only iia ^ (3. By assumption the 3-color instance (r^, r*^, r^, s^, s'^ , s^) 
is feasible, therefore the 2-color instance (r^, s^) is feasible as well — where yellow is renamed as 
black — which by Lemma [s] implies r^ ^ (s^)* and therefore also a ^ (3. This shows the first 
part of the lemma. 

Now assume that the instance has a solution, and a — f3. The n x {2n+ 2) grid is divided into 
3 parts (see figure [2]): into an n x n block (called RY-block), a nx 2 rectangle (called 2-column 
translator) and another nx n block (called GY-block). Again every block is sub-divided into an 
upper triangle, a diagonal and a lower triangle. 

Since a = (3, we have = (s^)*. So by Lemma [3] any solution must color in yellow the sJ 
first cells in every column j, and no other cell. In particular it means that the lower triangle of the 
RY-block must be red, the lower triangle of the GY-block must be green, and both upper triangles 
have to be yellow. 

Also on the first diagonal, the cell {i,i) has to be red if = 1 and yellow otherwise. On the 
second diagonal, the cell {n + 2 + must be yellow if = 1 and green otherwise. 
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Figure 2: The structure of the gadget (left) and a reahzation (right) for n = 7,k^3,u~2, v = 5 
and a = /3= (0,0,1,0,1,1,0). 



What can we say about the colors of the translator? If = a^, = 0, then the cells (n+1, u), {n+ 
2, u),{n + 1, v),{n + 2, v) have to be all red to satisfy the row projections. This contradicts the 
column projection s^_^_l = 1, and hence the instance is not feasible. 

Conversely, assume au + (Xy > 1. We will color the cells of the translator in a manner that 
respects the required projections. If i ^ {u,v} and ai = I — that is (z, i) is red — we color the 
cells (n + 1, i), (rt + 2, i) in green. If z ^ {u, v} and a.^ — 0, we color the cell (n + 1, i) in green and 
(n + 2, i) in red. 

Without loss of generality assume that a„ = 1. Hence {u, u) is red and we color (n + 1, u) in 
green and (n + 2, u) in red. We color {n + 1, v) in red. In addition we color (n + 2, v) in red if 

= and in green otherwise. It can be verified that the coloring defined above is a solution to 
the instance, which concludes the proof of the lemma. □ 



3 The reduction 

In this section we will construct a reduction from Vertex Cover to 3-COLOR TOMOGRAPHY. 
We basically use the same approach than in [2j, but with a different gadget. 

Vertex Cover is a well known intractable problem, indeed one of the first 21 problems shown 
to be NP-complete by Karp [B]. 

Vertex Cover Problem 

Input: a graph G — (V, E) and an integer k. 

Output: a set S C V oi size jS"! = /c such that for every {u,v) G E, u G S or v G S. 

Given an instance (G, k) of Vertex Cover, we construct an instance (r^, , , s^) of the 3- 
COLOR TOMOGRAPHY PROBLEM which is feasible if and only if the former instance has a solution. 
Without loss of generality we assume that fc < n — 2. 

Let be n = |T^|,m = \E\, and N — k{n ~ k){m — 1) + 1. We denote the m edges a,s E — 
{eg, ei, . . . , e,„_i}, and the n vertices as V — {1,2, ... , n}. We define an instance with N{n+1) + 1 
rows and N(n + 2) + n columns. 

For row p = 1, . . . , N{n + 1) + 1, let 

X- Lb -!)/("+ 1)J 
i = (p — 1) mod {n + 1). 

We think the set of rows as divided into N blocks of n+1 rows each, and a last block with a single 
row. We have x as the block index and i the row index relative to the block, with < i < n. Let 
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t = X mod m and consider the edge et — (it, v). We define the projections 



x{n + 2) 




if a; < and i = 
a X = N and i = 
if z G {u, v} 

.,n}\ {u,v} 

if x = and i = 
k if .X > and i = 
if z e {u, v} 

Hie {l,...,n}\ {u,v}. 




In the same manner, for column q = 1, . . . , N{n + 2) + n, let 

y^[{q-l)/{n + 2)\ 

j = ((q - 1) mod {n + 2)) + 1. 

The reason for defining j this way, is that if cell {p, q) is part of an RY-block or an GY-block, 
then will be the relative position inside the block with ranges I < i, j < n. Similarly as for 
the rows, we think the set of columns as divided into N blocks with n + 2 columns each and a 
last block with only n columns. Again, we have y as the block index, and j as the column index 
relative to a block with I < j < n + 2. For block 0<y<A^— Iwe define the red column 
projections as 



= (iV-2;-l)(n+l) + l + 



n-j + 1 if j e {!,..., n} 

1 if j = n + 1 

n-k+1 ifj = n + 2. 



For y 
are 



N we set — 0, for each j = 1, ... ,71. Similarly, for y = the green column projections 



and for I < y < N they are defined as 



n - 
n - 



5,^ = (y-l)(n+l) + l 



j a j e {!,..., n} 

n — I ifj = n+ l 
k -1 if j = n + 2. 



Clearly this a polynomial time reduction. It remains to show the following theorem. 

Theorem 1 The 3-color tomography instance is feasible if and only if the vertex cover instance 
is feasible. 

Proof: For one direction of the statement, assume that the vertex cover instance is feasible, and 
let b € Xn^k be the characteristic vector of a vertex cover of size k, i.e. 6^ = 1 if and only if i 
belongs to the vertex cover. 

We construct now a solution to the tomography instance. Consider the partitioning of the 
grid, as in figure |3] For convenience we refer to the source also as the 0-th row translator and to 
the sink as the {N + l)-th row translator. The j-th cell of the x-th row translator is defined as 
{x{n + 1) + 1, x{n + 2) + j). We color the R-frame in red and the G-frame in green. 

Let be any I < j < n. We color the j-th cell of the source in yellow if bj — 1 and in red 
otherwise. For a;=l,...,iV— Iwe color the j-cell of the x-th row translator in green if bj ~ 1 
and in red otherwise. In the sink we color the j-th cell in green if bj — 1 and in yellow otherwise. 

Now for block x — 0, . . . , N — 1, consider the instance to the gadget defined by a = /3 = 6, and 
u,v such that {u,v) = exraodm- By Lemma [i] it is feasible, since & is a vertex cover and hence 
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Figure 3: The general structure of our reduction 



bu + K > 1- Then we color the (n+l)x(2n+2) cells starting at ((a;-l)(n+2) + l, + 
exactly as in the solution to the gadget. It is straightforward to check that this grid satisfies the 
required projections, and therefore the tomography instance is feasible. 

For the converse, assume that the tomography instance has a solution. For every x = 1, . . . ,N 
we apply Lemma[T|for the red color and intervals / = [a;(n+l) + l, N{n+1) + 1] and J = [1, x{n+2)]. 
We deduce that in the solution the R-frame must be all red, and all GY-blocks (and also the G- 
frame) must be free of any red. Similarly, we show that the G-frame must be all green, and all 
RY-blocks must be free of any green. 

This implies that in the source, k cells are yellow, and n ~ k are red, in the row translators k 
cells are green and n ~ k red, and in the sink k cells are green and n — k yellow. We define the 
vectors 6°, 6^, . . . , 6^ e Xn.k, such that for all 1 < < n we have 

• = 1 iff the j-th cell in the source is yellow, 

• foj = 1 iff the j-th cell in the x-th row translator is green, for all 1 < x < iV. 

For X = 0, . . . ,N, consider the part P of the solution that is the intersection of rows [x{n + 
1) + 2, x{n + 1) + n + 1] and columns [x{n + 2) + 1, x{n + 2) + 2n + 2]. We number the rows of 
P from 1 to n and the columns from 1 to 2n + 2. Let {u, v) — mod m- By subtracting from the 
row projections the number of red and green cells in the frames, we deduce that row 1 < i < n 
in P contains i + I red cells and i green cells if z G {u, v} and i red cells and i + I green cells if 
i ^ {u,v}. 

We proceed similarly for the columns n + I and n + 2. By subtracting from the column 
projections the quantities that are in the frames, we deduce that column n + 1 of P contains one 
red cell, and n — 1 green cells, and column n + 2 contains n — k + 1 red cells and fc — 1 green cells. 

Column x{n + 2)+j for 1 < j < n contains n~j + l red cells that are not in the R-frame. Since 
GY-blocks are free of red, these cells must either be in the x-th row translator or in column j of P. 
Note that the j-cell of the x-th row translator is red iff 6J = 0. Therefore column j of P contains 
n — j + bj red cells and no green cell. Similarly column n + 2 + j of P contains n — j + 1 — b^^^ 
green cells and no red cell. 

This implies that P is the solution to the gadget defined by u, v, a, (3 with a = b^ and /? = b^^-^. 
Then by Lemma |4] we obtain that b^ ^ 6^+^ and in general 

6° ^ 6^ ^ . . . ^ 5^. 



6 



By the choice of N and Lemma |2] there exists an £ such that 

b' = 6^+1 = . . . = 6^+™. 

By Lemmajij we have b^+b^ > 1 for all {u, v) e {eg, ei+i mod m, ■ • ■ , ei+m~i mod m} = E. Therefore 
b^ encodes a vertex cover of size k, and this completes the proof. □ 

4 Related problems 

4.1 Edge-colored graphs with prescribed degrees 

We can reduce the 3-COLOR tomography problem to a similar graph problem. 

Finding edge-colored graphs with prescribed degrees Let be a set of two colors {i?, G}, 
and a vertex set V . We are given prescribed degrees dP ,dP : ^ N and have to find two disjoint 
edge sets E^,E'~^ C such that the graph G{V,E^ U E'~^) has the required degrees, i.e. for all 
veV 

d^{v) = \{u : {u,v) G d'^iv) = \{u : {u,v) € E^}\. 

Note that in contrast, finding an uncolored graph with given degree sequences can be solved 
in polynomial time, see for example |7]- 

Lemma 5 The problem of finding an edge-colored graph with prescribed degrees is NP-hard. 

Proof: We reduce from the 3-color tomography problem. Let (r^, , s^, s^) be an m x n-instance 
of the 3-color tomography problem. We set k = n + m, V = {I, . . . ,k}, and the following degrees, 
for 1 < i < n and 1 < j < m 

d^{i)^rf' + n-l d'^(z) = rf 

d^(n + j) = sf d'^{n + j) = sf + m-l. 



Now we show that the instance {r^,r'~^,s^,s'~^) is feasible if and only if the instance {d^,d'^) is 
feasible. For one direction, assume that there is a solution M to the 3-color tomography instance. 
We construct a solution E^, E'^ to the graph problem as follows. For any \ <i <n and I < j < m, 
a Mij = R, then {i,n + j) € E^ , if = G, then {i,n + j) e E'^ . Also for any 1 < i < i' < n, we 
have (i, i') € E^ and for any 1 < j < / < m, we have {n + j,n + j') S E*^ . Now clearly E^, E'^ 
satisfy the required degrees. 

For the converse, we define the quantity <i> = J^^^i d^i'i-) ~X]JLi d^{n + j). By assumption (jlj 
this value is n(n — 1). Since this value equals also 

|£:^n {!,..., n}2| - \E"n{n+l,...,n + m}\, 

there is a red edge between every pair of vertices (z, i') with I < i < i' < n, and no edge between 
every pair of vertices (n + j,n + j') with 1 < j < j' < m. Similarly we can show that there is a 
green edge between every pair of vertices (n 4- j, 4- j') with 1 < j < j' < m. 

Now let M be the m x n grid, with cell {i,j) colored in red if {i,n + j) S E^, and in green if 
{i, n + j) e E'^. By the degree requirements, M is a solution to the 3-color tomography instance. 

□ 
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4.2 Tiling Tomography 

Tiling tomography was introduced in ^ , and it consists of constructing a tiling that satisfies some 
given row and column projections for each type of tiles we admit. 

Formally a tile is a finite set T of cells of the grid N x N, that are 4-connected, in the sense that 
the graph G(T, E) is connected for E = {((i, j), (i',/)) : \i ~ ^'| + \j -f\ = 1}. By T + {i',f) = 
{{i + i' ,j + j') : & T} we denote a copy of T that is shifted i' units down and / units to the 
right. We say that a set of tiles is feasible if they do not intersect. In addition we say that it tiles 
the m X n grid if its (disjoint) union equals the set of all grid cells, and we refer it as a tiling. 

In the tiling tomography problem we are given a finite set of tiles T = {Ti , . . . , T^}, and vectors 
j.d g f^m^ g p^jn for 1 < (i < fc. The goal is to compute a matrix M e {0, 1, . . . , k} such that the 
set 

(J {Td + {i,j):l<i<n,l<j<m such that M,j ^ d, for 1 < d < A:} 

l<d<fe 

is a tiling of the m x n grid, with the projections 

rf = \{j : M,, = d}\ s'^ ^ \{i : Ah, = d}\. 

By width and height of a tile T we understand the size of the smallest intervals /, J such that 
T C / X J. This definition extends to set of tiles. A tile T is said to be rectangular if for every 
{i',j') such that {T,T + is feasible, we have that the width of {T,T + is at least 

twice the width of T or the height of the set is at least twice the height of T. 

It was conjectured in [3 , that for Ti being a single cell and T2 a non-rectangular tile, the 
{Ti, T2}-tiling tomography problem is NP-hard. This question is still open and intriguing. 

4.3 Rectangular tiles 

Consider two rectangular tiles, Ti being a pi x gi rectangle and T2 a, p2 x q2 rectangle, i.e. 
Tc — {0, . . . ,Pc — 1} X {0, . . . ,qc — 1}, for c G {1,2}. What can be said about the complexity of 
the {ri,T2}-tiling tomography problem? 

Ifgcd(pi,p2) = d> 1, then clearly any solution M S {0, 1, 2}™^" to a {Ti, T2}-tiling tomogra- 
phy instance (r^, s"^), must satisfy that if Mij ^ 0, then i mod d = 1. Therefore the {Ti,T2}-tiling 
tomography problem can be reduced to the {T[,T2}-tiVmg tomography problem, with T[ being 
a {pi/d) X qi rectangle, and a {p2/d) x q2 rectangle. We omit the formal reduction, which is 
straightforward. 

From now on suppose that gcd(pi,p2) = gcd{qi,q2) = 1. We distinguish the following cases, 
up to row-column symmetry. 

• If pi — p2 — 1, that is the tiles are two horizontal bars of length qi and q2, then the problem 
can be solved in polynomial time (Theorem |2]) . We use an idea already present in j4j , where 
it is proven for qi ^ I. 

• If pi = (72 = 1 and p2 = qi = 2, then the tiles are called dominoes, and again the problem 
can be solved in polynomial time, although with a more involved algorithm [5]. 

• If pi = 92 = 1, P2 > 2 and gi > 3 then the problem is open. The first author conjectures that 
the problem is NP-hard, while the other two conjecture that it could be solved in polynomial 
time with a similar approach as in [9]. 

• If pi, (?i > 2, then the problem is NP-hard (Theorem [3|. In [3j the special case pi = qi = 2, 
P2 — q2 — ^ was related to the 3-color tomography problem, and it is therefore also NP-hard. 
We generalize this reduction in section [4. 6| 

• If there is a third rectangular tile T3, then for the tile set {Ti,T2, T3} the problem is NP-hard, 
see section l4n 
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Figure 4: Transforming a solution M for c = 1. 



4.4 An algorithm for vertical bars 

Theorem 2 The tiling tomography problem can be solved in polynomial time for two rectangular 
tiles of dimensions pi x I and p2 x I. 

Proof: The algorithm is the simple greedy algorithm, as the one used in [4j. It iteratively 
stacks bars in the matrix. 

Formally the algorithm is defined like this. We construct a matrix A e {0, 1, 2}™^" with the 
required projections. Initially A is all 0. We maintain a vector v such that Vj is the minimal i such 
that Ai j =/: 0, and Vj = m + I if column j of A is all zero. Initially vj — m + I for all I < j < n. 
We also maintaing vectors f^,f'^,s^,s'^, which represent the remaining projections. Initially they 
equal the given projections of the instance. The vectors {v,f^,f^,s^,s'^) define a more general 
tiling problem, where in every column j, only the first Vj — 1 cells have to be tiled. 

The algorithm: Let i — maxvj. If i = 1 we are done, and return A, if all vectors 
f^,r'^,s^,s'^ are zero, and return "no solution" otherwise. 

If i > 1, let ii = « — pi and 12 — i—P2- If ^1-^ = = 0, abort and return "no solution" . 
Otherwise let c G {1,2} such that f^^ > 0. Let j be a column with Vj + \ = i that 
maximizes s^. Then drop the bar pcX 1 in column j, i.e. set Ai^j = c, and decrease 
ff^ and s^. Repeat the whole step. 

Clearly, if this algorithm produces a matrix, then it defines a valid tiling with the required 
projections. We have to show that if the instance has a solution, then the algorithm will actually 
find one. For this purpose, let be some step of algorithm such that the intermediate instance 
I := {v, r^, f^, s^, s^) is feasible. The initial step could be a candidate. Let M be a solution to it. 
Let i = maxvj. If i = 1, then f^, . . . , are all zero, since the instance is feasible. 

Let ii = i — pi and 12 = i — P2- We have that either j = 1 or Af^^.j = 2 for every column 
j satisfying Vj — i, since M is a valid tiling. Therefore some of r}-^,rf^ must be non zero. Let 
c, j be the values the algorithm chooses. Let T' be the instance obtained after the iteration of the 
algorithm, that is ff^, Sj are decreased by 1 and Vj by Pc- 

If Mi^ j = c, then M' which equals M except for Mi^ j = is a solution to T'. 

If Mi^j ^ c, then by the projections, there must be a another column k with Vk — i and 
Mi^^k = c. We will now transform M such that Mi^j — c. Then we are in the case above and 
done. 

By the choice of the algorithm we have sj: < Sj. By this inequality, there exists ia such that 
the total number of c's below the row is the same in both column j and column k. Take io 
being the largest one satisfying that. By the choice of io we have that Mi^j = c and Afi^ j c. 
Since M is a valid tiling, then the restriction to cells below iq in column k is also a tiling and then 
Migj 7^ 0. We conclude that between and i the number of I's and 2's in column j is the same 
as in column k. Then exchanging the parts of columns j and k in M between zq and i, does not 
change the projections of M, and we obtain the required property Mi^ j — c. 
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By the choice of the algorithm wc have < s^. Now we claim that there is a row ia such 
between io and i the number of I's and 2's in column j is the same as in column k. Indeed, 
consider the largest io such that the total number of c's below the row io is the same in both 
column j and column k. It must be that Mi^j = c and ^ ^ since M is a valid tiling. Then 
exchanging the parts of columns j and fc in M between io and i, does not change the projections 
of M, and we obtain the required property Mi^j = c. 

□ 

4.5 A general NP-hardness proof structure 

In the next section we will reduce the 3-color tomography problem to the tiling tomography 
problem for some fixed set of tiles T. The proof uses a particular structure that we explain now. 

Let (r^, r*^, r^, ,s^, .s*^, ,s^) be an instance to the 3-color tomography problem for an to x n 
grid. In the reduction we will choose constant size grid £ x k — that we call a block — and three 
T-tilings of it, that we denote M^, , . There will be two requirements: Let f^'"^, s'^''^ be the 
Td-projcctions of the tiling M'' for c e {R, G, Y} and d G {1,2}. 

The first requirement is that the vectors {f^'^'i, f^^'i, f^'^} are affine linear independent. 
The same requirement holds for the column projections {s^'^, s'^'^, s^'^}. This implies that every 
vector r spanned by f^^.i^ j^g^g ^ unique decomposition into r = hrt^ + ncf'^ + nyf^ 

for nR + no + ny = n. 

The reduction, consists of an mixnk grid, and the projections l<i<i, 'i-<j<k,l<x<m, 
1 < y < n, d e {1,2} 

'xi-e+i 

''yk-k+j 

The idea is that the m£ x nk is partitioned into mn blocks of dimension £ x k. The second 
requirement is that in every solution M to the tiling instance, all blocks of M, are either 
M^,M^ or blocks that have equivalent projections. 

Lemma 6 The instance to the T -tiling problem has a solution if and only if the instance to the 
3-color tomography problem has a solution. 

Proof: Let M G {i?, G,y}™^" be a solution to the 3-color tomography problem. We transform 
it into a matrix M G {0, 1, 2}™^^"*^ by replacing each cell (z, j) of M by the £ x k matrix M'^ for 
c = Mij. By construction, this is a solution to the tiling problem. 

For the converse, suppose that there is a solution M to the tiling problem. By the second 
requirement, every block of M can be associated to one of the colors {R,G,Y}. We construct 
a matrix M G {ii,G,F}'"^" such that M^y = c if the block {x,y) of M is M", or something 
projection equivalent. 

Fix some arbitrary 1 < a; < to. By the first requirement, the projections of the rows x£ — i -\- 
1,. . . ,x£ have a unique decomposition into njir^'^ + ncr'^-^ -\- nyr"^'^ with ur -\- no + ny = n. 
By the definitions of the projections riR = r^^riQ — r^,ny = , and then row x of M has the 
required projections. We proceed in the same manner for the columns and show that M is a 
solution to the 3-color tomography instance. □ 

4.6 An NP-hardness proof for two rectangular tiles 

Theorem 3 The tiling tomography problem is NP-hard for two rectangular tiles of dimensions 
Pi X qi and p2 x with gcd(pi,p2) = gcd(gi,g2) = 1 andpi,qi > 2. 
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Figure 5: The 3 valid block tilings. 



Proof: We apply Lemma [6] for £ = 2pip2 and k = 2qiq2. The 3 tilings of the £ x k grid are 
depicted in figure |5] and defined formally as follows. The rows / = {!,...,£} and the columns 
J — {1, . . . ,k} are partitioned into sets /i, /2, /s, I4 and Ji, J2, J3, J4 defined as 



h = {l,...,P2} 

h = {P2 + l,---,PlP2} 

-^3 {PlP2 + !,... ,PlP2 +P2} 

^4 = {PlP2 + P2 + 1, • • ■ , 2piP2} 



Jl = {1, . . . ,(?2} 

J2 = {Q2 + 1, • • ■ ,'?1'?2} 

-^3 = {qiq2 + 1, . • . ,<7i'72 + 92} 

J4 = {<7i(72 + 92 + 1, ■ ■ . ,291 92}- 



Then AI^ is defined as the block tiling that covers (/i U I4) x ( J3 U J4) with r2 and the rest with 
Ti, M*^ is defined as the block tiling that covers (Is U J4) x (Ji U J4) with T2 and the rest with 
Ti, while is defined as a tiling using only Ti. These tilings are uniquely defined. Clearly the 
row Ti-projections of the 3 tilings are affine linear independent, so the first requirement of the 
construction is satisfied. 

The second requirement follows from a sequence of observations. Let M be the solution to the 
tiling instance, obtained by reduction from a 3-color instance {r^ , r^, r^, , s^, s^). 

First note that in the tilings M^, M*^, , every tile is completely contained in the £xk block. 
Therefore the tiling instance has zero projections for Ti at rows x with (x — 1) mod £ > £ — pi + 2. 
A similar observation holds for tile T2 and for the column projections. As a result in M every tile 
is completely contained in some £x k block, and in other words every block of M is {Ti, T2}-tiled. 

What can we say about the possible tilings? Again note that in the tilings M^,M'^,M^, 
every row in I2 is completely covered by Ti-tiles. Therefore by the projections, this holds also for 
every block in M. The same observation can be done about columns in J2- 

Note that if api+bp2 = 2pip2, then (a, b) £ {(0, 2pi), {2p2, 0), {p2,Pi)}- This is simply because 
by gcd(pi,p2) = 1, in any solution to api — p2{2pi — h), a must be a multiple of p2. Together with 
the previous observation, this implies that every column of a block is either covered completely by 
Ti-tiles or covered half by Ti-tiles and half by r2-tiles. The same observation holds for the rows. 

The trickiest observation of this proof is that in every block of M, the region Ii x Ji is covered 
by Ti. For a proof by contradiction, suppose it is covered by T2, in fact by a single tile T2 since 
|/i X Ji| = IT2I. But since /2 x J is covered with T2, and by gcd ((71,^2) = 1, it must be that the 
cell (p2 + 1,<72 + 1) is covered by a tile T2 + {p2 + 1, j) for some column j < 52 ■ By the same 
argument, the cell {p2 + 1, ^2 + 1) is also covered by a tile T2 + (^2 + 1, i) for some row i < p2- 
Therefore these two tiles overlap in {p2 + 1, 52 + 1), which contradicts that AI is a (valid) tiling. 

Now fix a block of M. If row 1 is partly covered by T2, then T2— tiles must cover the half 
columns in J. Hence in the row 1 they cover exactly the columns in J3 U J4. The same argument 
shows that every column j G J3 U J4 is then half covered by r2— tiles. Previous observation 
state that I2 x {j} is covered by Ti. But the length of I2 is a not a multiple of pi. Therefore 
{PiP2 + 1, j) must then also be covered by Ti and hence {I2 U 13) U {j} is covered by Ti— tiles. 
Therefore (Ji U I4) x {j} is covered by T2. The choice of j was arbitrary, and therefore the 
block-tiling is exactly . 
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Similarly we deduce that if column 1 is covered partly by T2, then the block-tiling is exactly 
M'^. Now if row 1 and column 1 are completely covered by Ti, then (/i U x J and / x ( U J2) 
are completely covered by T2— tiles. As a result the block-tiling only contains in (/3U/4) x (J3U J4) 
either Ti— tiles or tiles, that correspond with the tiling and another we call the bad tiling, 
respectively. 

We will show that no bad tiling appears in M. Let Nu be the number of blocks in M that are 
M^. Similarly, let Nb be number of bad block-tilings in M. Note that the row projection of a 
bad tiling equal the row projections of and that the column projections equal the projections 
of M^. Therefore by the projections we have the equalities 

i j 

Since by assumption J^i^f — J2j ^f' have Nb = 0. This shows the second requirement of our 
construction, and by Lemma |6] completes the proof. □ 

4.7 An NP-hardness proof for three rectangular tiles 

Theorem 4 The tiling tomography problem is NP-hard for any 3 rectangular tiles. 

Proo/.' [sketch] Let pi x qi, p2 x q2 and ps x the respective dimensions of 3 tiles Ti,T2,T^. 



The idea of the construction is that we apply the general proof scheme from section 4.5 with 3 
tilings M^,M'^,M^, such that contains tile Ti in position (0,0), M° contains T2 and 
contains in position (0, 0). Moreover each of the 3 tiling minimizes lexicographically rii, n2, n^, 
where Uc is the number of tiles in the tiling. 

Formally, let ii be the smallest number i > with i modpi = and either i modp2 = or 
i mod p3 = 0. Let i2 be the smallest number i > ii with (i — ii) mod P2 = and i mod p^ = 0. 
Let is be the smallest number i > with i modp2 = and i modpa = 0. We define numbers 
ji,j2,i3 in exactly the same manner with 91,(72,93 playing the same role as pi,p2,P3- 

We apply Lemma for k — max{z2,J3} and £ = max{j2, J3}. The 3 tilings of the £ x k grid 
are depicted in figure (GT and defined formally as follows. In this section we assume for convenience 
that the rows and column indices relative to a block start at instead of 1 . 

In M^, the [l,ji] x [l,ji] subsquare is completely tiled with Ti. Then the region [1,12] x 
[lij2] — [Ij^i] X [l,ii] is completely tiled with T2 and the remaining part with T^. Note that by 
the choice of k, no column is intersects a tile Ii,T2 and T3. For example column intersects Ti 
tiles from row to ii — 1, and then either tiles T3 from row ii to k (if ii mod p3 = 0) or tiles T3 
from row ii to k (if ii mod pa 0). The same holds for rows. 

In M'~^ , the [1,«3] x [1,^3] subsquare is completely tiled with T2, and the remaining part with 
T3. In M^, the whole [l,fc] x [1,£] block is tiled with T3. 

Clearly the projections of these tiles satisfy the first requirement for the general proof structure. 
Now there are different interesting observations to make. In the tilings above, every tile T^, + {i,j) 
satisfies i mod Pc = 1 and j mod gc = 1 • This means that the projections of tile Tc are zero for 
any row i ^ 1 (mod Pc) or column j =/= I (mod qc), and in any solution M to the tiling instance 
resulting from the reduction, the property must hold for all blocks as well. This observation is 
crucial for the proof. 

In particular it implies the following fact (*). Fix some solution M to the tiling instance 
resulting from the reduction. Consider a block in M. If there is some vertical separation between 
two types of tiles, in the sense that cell (z — is covered by some tile Ta and (z, j) by some tile 
Th with a b, then we must have i mod Pa — ^ and i mod pb — I. 

We use this observation to show that the construction satisfies also the second requirement. 
Fix some solution M to the tiling instance resulting from the reduction. We distinguish 3 types 
of blocks: (1) Blocks that contain a tile Ti, (2) blocks that do not contain any tile Ti, but contain 
a tile T2, and (3) blocks that are completely tiled with T3. So blocks of the third type are exactly 
, and we have to show that blocks of the first type are exactly and blocks of the second 
type are M*^. 
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Pi = 2, gi = 5 p2 = 3, (J2 = 3 p3 = 5, 93 = 2 



Ti 



n 

T2 



T3 



ji = 15 j2 = 18 ia = 6 



ii — 6 



12 — 15 



is = 15 



Figure 6: The 3 valid block tilings. 



Let ri/j be the number of blocks of the first type, and consider one of them. Then by the 
projections all tiles Ti must be contained in the region x and by the observation (*) 

above, the whole region must be tiled with Ti. Now if i2 > H, then the region [1,«2] x [1,^/2] ^ 
X cannot contain tiles Ti nor T3 and must be completely tiled with T2. Later we will 

show that the remaining part of the block is tiled with T3. 

Note that nn is also the total number of red projections in the original 3-color tomography 
instance, so if 13 < «2, then all r2-tiles that have to be placed in a block-row i > max{i2, is}, are 
placed in a type 1 block. The same observation can be made for columns. Therefore all tiles T2 in 
a type 2 block, must be placed at positions of the form S [1, 13] x [1, J3]. By the observation 
(*), the whole region [1,13] x [l,j3] is completely tiled with T2. By the observation above, the 
remaining part can only be tiled with T3, which shows that the type 2 blocks are exactly M'~^. 

Let no be the number of type 2 blocks in M. It is also the total number of green projections 
in the original 3-color tomography instance. Therefore by the projections, every tile T2 + 
with i > max{i2,i3} or j > max{j2, J3} must be contained in a type 2 block. This shows that the 
remaining part [1, k] x — [1,«2] x [1,^2] of a type 1 block, contains only r3-tiles. This shows 
that the type 1 blocks are exactly M^. 

Therefore the construction satisfies the second requirement for Lemma [6] and we are done. 

□ 
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